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The lunar microcrater phenomenology is described. The morphology of the lunar craters 
is in almost all aspects simulated in laboratory experiments in the diameter range from 
less than 1 p. to several millimeters and up to 60 km/s impact velocity. An empirically de- 
rived formula is given for the conversion of crater diameters into projectile diameters and 
masses for given impact velocities and projectile and target densities. 

The production size frequency distribution for lunar craters in the crater size range from 
1 n to several millimeters in diameter is derived from various microcrater measure- 
ments within a factor of up to 5. 

Particle track exposure age measurements for a variety of lunar samples have been 
performed. They allow the conversion of the lunar crater size frequency production 
distributions into particle fluxes. These fluxes derived from lunar samples are consistent 
with satellite in-situ measurements. For example, for masses m > (1 ± 0.4) x 10' 2 g, the 
latest HEOS satellite flux of <j> = (7 ± 2) X 10~ s nr 2 s' 1 agrees with <t> = (1 — 7) X 10~ 5 m 2 
s" 1 derived from lunar data. 

The development of crater populations on lunar rocks under self-destruction by subse- 
quent meteoroid impacts and crater overlap is discussed and theoretically described. 
Erosion rates on lunar rocks on the order of several millimeters per 10 7 yr are calculated. 

Chemical investigations of the glass linings of lunar craters yield clear evidence of 
admixture of projectile material only in one case, where the remnants of an iron-nickel 
micrometeorite have been identified. 


Microcraters created as a result of hyper- 
velocity impacts of interplanetary particles 
are found on most lunar rock surfaces. In this 
respect the microcraters can be considered 
as micrometeoroid detectors exposed for geo- 
logical periods of time. 

It is the purpose of this presentation to 
outline the principal phenomenology of such 
craters and to give some synthesis of the in- 
vestigations of a working group at the Max 
Planck Institut fiir Kernphysik, Heidelberg, 
which is concerned with crater statistics, lab- 
oratory cratering experiments, solar flare 
track exposure ages, calculations of inter- 


planetary dust fluxes, erosion processes, and 
chemistry of craters (see also ref. 1) . 

Crater Phenomenology 

Microcraters on silicate target materials 
consist of a central glass-lined depression 
(central pit; center of impact) embedded in 
a so-called spallation zone where target mate- 
rial is spalled off under the influence of a 
shock wave propagating in the course of the 
impact process. On glass targets the spalla- 
tion zone tends to be circular and concentric 


585 


586 


COSMOCHEMISTRY OF THE MOON AND PLANETS 



Figure 1 . — Scanning electron microscope pictures of 
microcraters on lunar sample 12063,106. (a) and 
(b) Craters in crystalline material, (c) and (d) 
Craters in the glass linings of larger craters. 


to the central depression. Spallation zones on 
crystalline material are mostly irregular in 
shape and may follow crystal structures or 
grain boundaries. Occasionally the material 
affected by the shock process is not com- 
pletely thrown out of the spallation zone, but 
parts of it rest there slightly lifted above the 
undisturbed target surface (refs. 2 and 3). 
Examples of craters in glassy and crystal- 
line lunar material are shown in figure 1. 

In most cases larger craters ( > 50-/im cen- 
tral pit diameter) show a circular zone of 
higher albedo situated between the central 
pit and the spallation zone. This “halo zone” 
consists of finely fragmented material under- 
lying the central pit. A crater with a prom- 
inent halo zone is displayed in figure 2. 

In some cases the central pit is found to sit 
on some kind of a pedestal raised above the 
spallation level (“stylus pit”), or it may be 
not present at all (“pitless craters”) (refs. 2, 
4, and 5). An example of a stylus pit crater 
is given in figure 3, together with its labora- 
tory-produced counterpart. 

Microcraters on metal targets display only 


a central pit, surrounded by a raised lip of 
“impact-fluidized” material (refs. 6 and 7) . A 
similar behavior is found for very small 
( <. l-/xm diameter) craters on silicate mate- 
rial. The spallation zone is not present, and a 
depression looking very much like craters in 
metal targets remains. Often these craters 
are elongated, reflecting the obliquity of im- 
pact or irregularity of the impacting particle 
shape (ref. 8). 

The morphology of microcraters on lunar 
samples can be simulated by laboratory ex- 
periments in almost all aspects (ref. 4). Fig- 
ures 4 and 5 give a selection of lunar- and 
laboratory-produced craters with correspond- 
ing morphologies. 



Figure 2 . — Crater with prominent halo zone of higher 
albedo on lunar sample 60015,35 (central pit di- 
ameter 300 pm). 



Figure 3. — “Stylus-Pit" craters, (a) Lunar sample 
1202U,8,1. (b) Laboratory produced crater, 3 km/ 
s < v < 5 km/s, Af — » quartz glass. 
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Figure 4 . — Types of micron- and submicron-sized 
craters found on lunar samples, in comparison to 
craters produced in the laboratory. Lunar sample 
number: (A) 12063,106a Top (cryst. rock); (B) 
14257 (coarse fines) ; (C) 14257 (coarse fines) ; and 
(D) 12063,106a Top (cryst. rock). Impact con- 
ditions: (a) Fe— > Norite, v — 7.6 km/s; (b) Al -* 
Quartz glass, v = 5 km/s; (c) Fe — * Duran glass, 
v = 50 km/s; (d) Fe —* Duran glass, v = 30 km/s, 
impact angle 45° . 



50 urn ' a»70° 1 100 urn ' a»70° 1 3 pm 1 a«60° 


a b 

Figure 5. — (a, first two rows) Microcraters on Fe- 
Ni sample 60502,17 in comparison with (b, third 
row) laboratory -produced craters on metal targets 
(SEM-photographs, a = observation angle). 


Laboratory Simulation Experiments 

Using an electrostatic dust accelerator, 
projectiles in the mass range between 1(H 0 
and 10“ 16 g were accelerated to velocities be- 
tween 1 and 60 km/s (ref. 4). Using a lith- 
ium plasma gun, projectile masses between 
10~ 5 and 10 _8 g were accelerated to velocities 
between 2.5 and 18 km/s (ref. 9). Light gas 
gun projectiles of 10~-g were accelerated to 
velocities between 1 and 4.5 km/s (ref. 10). 
Projectile materials were alumuinum, carbon, 
iron, glass, and stainless steel, and target 
materials were glass, norite, stainless steel, 
aluminium, copper, and gold. 

For both the dust accelerator and the light 
gas gun, it results for normal impacts occur- 
ring throughout the velocity range that the 
ratio of pit diameter to projectile diameter, 
Dp/d, is proportional to the two-thirds 
power of the impact velocity. Thus, D P /d 
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Figure 6 . — Dependence of the central pit diameter, 
Dp, on projectile mass, m, or diameter, d, for dif- 
ferent velocities. Data are for micron-sized alumi- 
nium projectiles impacting stainless steel, norite, 
and glass. 

v 2/3 . In figure 6, the dependence of D P on d 
is given for different impact velocities. 

Plasma gun experiments, however, pro- 
duced mainly pitless craters. The ratio of the 
spallation zone diameter to projectile diame- 
ter, DJd, is proportional to the two-thirds 
power of the impact velocity; thus, D,/d 
~ v 2/3 . 

Though over-limited mass ranges, the ex- 
perimental results can in each case be ap- 
proximated by the two-thirds power laws, 
and the ratio D P /d really slightly increases 
with increasing mass. From experiments 
showing pits, D P /d increases a factor of 2 as 
projectile mass increases from 10 _12 g to 

10- 2 g. 

The ratio of the spallation zone diameter 
to projectile diameter, DJd, increases by 
a factor of 2 as projectile mass increases 
from 10 -12 to 10 _5 g. The ratio of the ki- 


netic energy of the projectile to the volume 
of the crater is E/ r — 3.6 X 10 11 erg/cm 3 for 
a projectile mass of 10~ 12 g, and E/ r = 4X10 10 
erg/cm 3 for a projectile mass of 10~ 2 g. 

Using a formula from Gault (ref. 11) 
given for micrometer-to-centimeter-sized 
projectiles, these results can be represented 
approximately for normal impact over the 
whole experimental range from micron-sized 
to millimeter projectiles by slightly changing 
the constants 

D 8 = KH 077 X P r 1/2 X p,, 1 ' 3 X E 1/2 - 65 

where p t and p p are target and projectile 
densities, respectively. 

Crater Statistics 

Cratered rock surfaces can display three 
categories of crater size frequency distribu- 
tions. 

As long as every impact crater on a sur- 
face exposed to the micrometeoroid bombard- 
ment is preserved, it is considered to be in 
the production state. The craters lie far apart 
from each other and do not superpose one 
another. Figure 7 shows such a production 
crater population on a glass-coated lunar 
rock. 

As time goes on, the crater frequency in- 
creases and crater superposition occurs; i.e., 



Figure 7 . — Production crater population on lunar 
sample 60015,35. The light bar on top corresponds 
to 1 mm. 




Figure 8 . — Closeup photograph of rock 66075,0, con- 
sidered to represent an equilibrium population. 
Width of frame approximately 1 cm). For com- 
parison, a sketch depicting actual central pit lo- 
cations and size (black dots) is added, and the 
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Figure 9 . — Cumulative frequencies of microcraters 
obtained on a variety of glass surfaces represent- 
ing production populations for microcraters with 


actual, as well as eroded, spall zones (dashed cir- 
cles) are indicated using the measured D S /D P 
value. Note the close proximity of many events, 
demonstrating that the effective destruction diame- 
ter operating in this surface is D s . Note also the 
relative paucity of small events within the ( some- 
what younger) spall zone of the big crater in the 
center. 


pit diameters >25 pm. The consistent data for 
1205U and 60015 are considered to be the most re- 
liable because of the large number of craters 
counted in both cases. Note the relative steepening 
of the production slope at about 200-pm pit diame- 
ter, resulting in a production slope steeper than 
—2. The curves are normalized at a pit diameter 
of 100 pm. 






UJO) N 


590 


COSMOCHEMISTRY OF THE MOON AND PLANETS 


diameter. The slopes for D > 250 pm are 
grouping around —3 and steeper. The differ- 
ences in slope are attributed to different ex- 
posure geometries of the individual lunar 
samples (ref. 12) and observational and sta- 
tistical uncertainties. The flattening of the 
curves for D < 250 /un is real, though some- 
what exaggerated by observational losses in 
some cases and influenced by geometry ef- 
fects. It means a depletion of craters of this 
size. 

Figure 10 shows cumulative production 
state crater number densities in the micron 
and submicron size range, where only scan- 



Figrure 10 . — Cumulative frequencies of micron-sized 
craters (production populations ) as a function of 
the central pit diameter. 


ning electron microscope investigations are 
feasible. These curves are not normalized and 
thus reflect differences in exposure age (ref. 
6 ). 

Measurements by different authors provide 
sufficient data to allow one to determine the 
production size distribution of microcraters 
on lunar rocks from less than l/i to almost 
1 cm in diameter. To derive the size distri- 
bution, we consider the following (1) Size 
distributions can usually be measured only in 
certain size ranges that depend on the size 
of a lunar rock, its exposure age, and the 
smoothness of its surface. (2) The absolute 
frequencies of craters in a certain size range 
vary with the exposure age of a rock. (3) In 
order to obtain the size distribution over the 
whole range of crater sizes observed on lunar 
rocks, normalization of microcrater frequen- 
cies derived from different measurements on 
different rocks is necessary. (4) The normal- 
ization range has to be chosen according to 
the quality of the data and the amount of 
overlap of the diameter ranges. 

Figure 11 gives a compilation of currently 
available measurements of microcraters on 
lunar rocks, reduced in the above discussed 
way (ref. 17). The size distribution is de- 
termined within a factor up to 5. The slope 
is approximately equal to —2 in the micron- 
to submicron-size range ; i.e., the distribution 
is steepening for D < 1 /un (see also ref. 18) . 
It is close to — 1 in the range 1/im < D < 100 
/im. In the range D > 300 /un, the slope has 
been determined to have the value —3.6 + 0.3. 
Thus, the distribution steepens significantly 
for the millimeter-sized craters. 

Another factor influencing crater statis- 
tics is the number of pitless craters that are 
mostly not accounted for. There is some evi- 
dence that the total number of pitless craters 
observed on cratered surfaces depends on 
mechanical surface properties, i.e., on the ab- 
solute crater frequency which controls the 
state of erosion. Although the statistics are 
relatively poor — eight data points for brec- 
cias and four data points for crystalline ma- 
terial — there is an indication that the num- 
ber of pitless craters decreases with increas- 
ing areal crater density, an observation 
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Figure 11 . — Compilation of microcrater measurements of different authors. The frequencies are normalized 
at the ranges indicated (vertical shift of the distributions). 


which can be explained in different ways. 
However, there is no unambiguous explana- 
tion yet available, so we prefer not to go 
into further details (see ref. 5). 

Superposition of craters changes the slope 
of the crater distribution. For any production 
distribution slope steeper than —2, it has 
been shown that the slope expected for the 
equilibrium crater population is —2 (ref. 
19). In figure 12 the development of a crater 


distribution with time is schematically 
shown. The longer the exposure time, t, the 
larger the craters affected by superposition. 
Thus, the production size distribution bends 
over to the equilibrium size distribution with 
the —2 slope for these sizes. With these as- 
sumptions, the tail of the distribution is al- 
ways in the production state. It reflects the 
differences in exposure time. An alternative 
view also exists (ref. 20). Because rocks 



592 


COSMOCHEMISTRY O* THE MOON AND PLANETS 



Figure 12 . — Schematic display of the development of 
a crater population. Larger and larger craters are 
destroyed by superposition as time, t, goes on. The 
tail of the distribution remains in the production 
state, whereas small craters are in the equilibrium 
state. 

destroyed early in their exposure lifetime are 
not collected by the astronauts, we can ex- 
pect to recover a disproportionately large 
number of rocks which have survived longer 
than the expected mean lifetime. Thus, most 
rocks collected may not have received a rep- 
resentative share of larger craters and have 
enjoyed a statistical selection effect. Those 
large craters not accounted for would cause 
an artificial steepening of the cumulative 
size distribution curve. Therefore, lunar rock 
surfaces that have already reached an equi- 
librium state for small craters do not neces- 
sarily yield valid production distributions for 
larger craters. Shown in figure 13 is a hy- 
pothetical log-log cumulative size distribution 
diagram having a reference equilibrium line 
with a slope of —2. Additional lines are shown 
that indicate how the observed distribution 
would appear for different numbers of unac- 
counted large craters due to the selection ef- 
fect. 


Figure 14 shows the behavior of real popu- 
lations on lunar rocks (ref. 17). The picture 
is similar to the schematic treatment in fig- 
ure 12. The difference in the frequencies at 
the tails of the distributions is greater than 
10. The frequencies tend to become equal for 
smaller sizes where the slope is approximately 
equal to —2. The difference between the high- 
est and the lowest frequencies is less than 
a factor of 3. Ideally, the frequencies should 
be equal. This spread in the equilibrium fre- 
quencies is explained by different destruc- 
tion conditions because of different ratios of 
spall sizes to pit sizes on different rocks 
(ref. 12), irregular spallation, and unob- 
served pitless craters. Additionally, the low- 
frequency populations might be still in the 
transition state. 

Quantitatively, the equilibrium frequency 
N t: is writted as N s = A X D -, where the 



Figure 13 . — Diagram illustrating the steepening of 
an equilibrium size distribution curve due to selec- 
tion effects. 
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crater diameter D can be either the spalla- 
tion diameter D, or the pit diameter D P . 
T he areal density coefficient, A, is expressed 
in the following way: 


A 


2 X (« -2) 

7T 



In general, however, crystalline rocks seem 
to have lower maximum densities than brec- 
cias. 


Exposure Ages 


where — a is the slope of the production pop- 
ulation in the log N versus log D plot. The 
Dp/D, ratio is generally smaller on crys- 
talline rocks than on breccias. Thus, one has 
to expect a corresponding variation in ab- 
solute equilibrium crater frequencies as is 
demonstrated in figure 14. 

The highest frequencies we found on lunar 
rocks are limited by the above-discussed pro- 
cesses. Maximum crater number densities 
found per square centimeter are 100-150 for 
D P > 100 /un; 3000 for D P > 1 /im; 10* for 
Dp > 0.1 n m (refs. 6 and 12). 

As discussed above, the density of craters 
can be expressed independent of size if the 
crater distribution follows the law N = 
A X D~ 2 . This is valid for populations in 
equilibrium and for the range D P < 250 /un 
as a tangential relationship for any lunar 
microcrater population because the slope flat- 
tens for smaller sizes. 

Maximum measured areal density coeffi- 
cients for a number of different rocks with 
different physical properties are plotted in 
figure 15, together with the limiting fre- 
quency value of Morrison et al. (ref. 21). 
Gault’s (ref. 22) saturation percent values 
are given for comparison with the areal den- 
sity coefficient values (ref. 12). We find, es- 
pecially for the Apollo 16 rocks, a wide range 
of crater densities. Low maximum values for 
“glass” surfaces occur because we measured 
essentially populations in the production 
state. Other populations with similar values 
may be considered also in production, pro- 
vided DJDp is the same, and excessive num- 
bers of pitless craters do not exist. 

Regarding the higher areal density values, 
we have concluded that these populations are 
in or approaching the equilibrium state. We 
find no simple relationship between areal 
density values for populations close to equi- 
librium and physical properties of the rocks. 


In order to convert crater populations on a 
given sample into the interplanetary dust 
flux, the sample’s residence time at the lunar 
surface must be known. Therefore, solar 
flare exposure ages were determined for sam- 
ples 15076,31; 15015,24; 15927,3; 15301,79; 
15205,51,1; 60502,17; 60015,35; and a glass- 
iined crater pit found in the Luna 16 core 
(C 19.118 ; 20-22 cm) . For track observation, 
both optical and scanning electron micros- 
copy were used. Detailed procedures are de- 
scribed elsewhere (ref. 23). 

In all but samples 15076 and 60015 the 
depth dependence of solar flare tracks was 
studied in the glass phase. In lunar glasses 
the solar flare tracks are annealed to a con- 
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siderable extent (refs. 4 and 24). Accord- 
ingly, the track sizes are diminished and the 
track densities reduced (ref. 24). However, 
this very fact allows accurate track counting 
in a sample’s near-surface, high-track den- 
sity region and evaluation of the depth de- 
pendence of solar flare tracks. Thereafter, 
the track densities were corrected according 
to the respective size shrinkage of the track 
etch pits (ref. 24) and/or normalized to the 
track densities of pyroxenes that were 
trapped within the glass phases. 

In all samples studied (except sample 
60502, probably due to a 100-/<.m Ni-Fe 
cover) the track densities P decrease with 
depth R from the surface according to P = 
const X R~ a . The «-values range between 2.5 
and 2.8 (ref. 23), agreeing well with the 
track record of the uneroded Surveyor III 
glass lens (refs. 25, 26, and 27). We there- 
fore conclude that erosion affected our sam- 
ples only to a minor degree. Furthermore, 
the surfaces were irradiated and cratered 
having a solid angle of exposure of essen- 
tially 27 r. The solar flare exposure ages are 
compiled in table 1. For their calculation, 
track production rates based on the revised 
Surveyor III track record were used (ref. 
23). (In an earlier paper, ref. 6, the track 
production rates used were different and the 
ages therefore lower.) The Surveyor III 
spacecraft was on the Moon during the pe- 
riod of maximum solar activity and for a 
fraction of the present solar cycle only. On 
the other hand, the flux data have not to be 
representative for the long time averaged 
solar activity. This uncertainty in applying 


the Surveyor III record to our lunar samples 
results in an age uncertainty of about a fac- 
tor of 2. 

Among charged particle track exposure 
ages of lunar rock surfaces determined by 
different groups, there exist some apparent 
systematic discrepancies up to a factor of 3. 
This might be due to an insufficient calibra- 
tion and standardization of the particle track 
method. For example, there exists a differ- 
ence in the etchable ranges of tracks in the 
various detectors used (glasses, feldspars, or 
pyroxenes), differences in the track etching 
conditions applied (etchant, temperature, 
and etching time) , and differences in the 
track observation technique (optical, scan- 
ning electron, or transmission electron mi- 
croscopy) . 

Particle track exposure ages have been de- 
termined by different authors for a variety 
of Apollo 12, 14, 15, and 16 rocks, where 
mostly millimeter-sized craters have been 
measured. As seen before, the size distribu- 
tion in this range follows the law 
N Dp- 3 57 on the average, or in terms of 
spall diameters, N = Q X D,- 3 87 . Normalized 
to crystalline rock, we obtain N Borm = <2 norm X 
Df 3 57 , where Q norm = q 3 - 57 X Q. The factor q 
is the ratio of DJD P for crystalline rock and 
D s /Dp for the special rock whose crater popu- 
lation is measured: q — {?* ^ crys ^ . It is as- 

D„ (meas.) 

sumed that for the same impact on different 
rocks with different composition and texture 
(e.g., breccias compared to crystalline rocks), 
the central pit with diameter D P has the same 
size, but the ratio of spall diameter to pit 


Table 1 . — Exposure Ages Determined by Means of Solar Flare Tracks 


Sample number 

15015, 

24 c.> 

Luna 

16"» 

15927,3 


HI 


60502,17'” 

60015,35 

Solar flare exposure 
age 

40 

5.7 X 10 2 

4.8 X 10’ 

7.9 x 10‘ 

9.6 X 10* 

2.6 X 10= 

1.1 X 10‘ 

1.5 X 10= 


Notes: (1) No microcraters have been found on these samples. 

(2) Considering the higher stopping power of iron, we corrected our data to an equivalent shielding 
of stony material. 
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diameter, D,/D P , is different. The time inte- 
gral of the cumulative crater production rate 
gives N = <f>(D) X t A = const X t A 
X D,- 3 57 (time constant influx assumed). 
Thus it results after normalization 

(a) for D, = 1 mm 

N norm ^1 X £4 

or 

(b) Q norm = c t X t A , independent of D,. 
Particle track exposure ages of dated rocks 

are plotted against the respective crater fre- 
quency measurements in terms of N„ orm (D a — 
1 mm) and Q norm in figure 16. The data can be 
fitted by the expression 

t A = 7.3 X 10 7 X Qnorm 
respectively, 

t A = 2.10 4 • N n orm (D» = 1 Him) 

which would correspond to a time constant in- 
flux between roughly 3 X 10 5 and 3 X 10® yr 
ago. These formulae can be used to determine 
exposure ages by crater statistics that corre- 
spond to particle track ages. 

A different view is presented by Hartung 
et al. (ref. 20). The correlation shown in 
figure 16 may be independent of time be- 
cause both solar flare tracks and microcra- 
ters occur at the surfaces of lunar rocks. 
Thus, rocks having physical properties per- 
mitting a higher erosion rate would have 
lower track and crater densities at the sur- 
face. The correlation observed then reflects 
different physical properties of rocks rather 
than different exposure ages. 

Meteoroid Fluxes 

Using the following formula it is possible 
to calculate the meteoroid flux </> for a pro- 
duction crater population with known expo- 
sure time. 


where </> -- meteroid flux per m 2 s 

N = cumulative crater number per m 2 
t = exposure time 


Knowing the impact velocity, the conver- 
sion of crater diameters into particle diam- 
eters is accomplished with the help of the 
simulation results. However, the results ob- 
tained in this way are for an isotropic 
spatial distribution. The results of recent 
space probe and satellite measurements show, 
however, a rather strong anisotropy in spa- 
tial distribution. 

The Pioneer 8/9 dust experiment (ref. 35) 
shows that a sensor orbiting the Sun at 1 AU 
and spinning in the ecliptical plane is hit 
by two major dust streams, one from the 
apex direction (direction of Earth’s motion 
around the Sun) and one from the direction 
of the Sun. Contrary to the apex particles, 
the Sun particles are the most frequent ones 
in terms of event rates. Although it was not 
possible to measure particle mass and impact 
velocity separately, there is strong evidence 
that these particles are extremely small 
(< 1 /im in diameter), and the impact ve- 
locity is higher compared to the apex parti- 
cles. The impact velocity is considered to be 
higher than 50 km/s and may even exceed 
100 km/s (ref. 35), depending on particle 
diameter. Zook and Berg (ref. 36) propose 
that these submicron-sized particles are being 
produced from normal zodiacal dust parti- 
cles. These zodiacal dust particles, following 
the Poynting-Robertson effect (ref. 37), are 
continuously losing energy and thus spiral- 
ing into the Sun. In the neighborhood of the 
Sun the particles partly melt and/or are 
fragmented by collisions. As a result of these 
processes, submicron-sized fragments are ac- 
celerated by the radiation pressure and blown 
out of the solar system on hyperbolic orbits. 
Hemenway et al. (ref. 38) claim, however, 
that submicron-sized particles may come off 
the surface of the Sun and may be blown 
out of the solar system. 

The HEOS 2 dust experiment (refs. 39 
and 40) is measuring the dust population in 
apex, anti-apex, ecliptic north, and ecliptic 
south directions. It has shown that the apex 
flux is the most dominant flux. The fluxes in 
the other directions are approximately 2 or- 
ders of magnitude lower compared to the 
apex flux. The impact velocities of the apex 
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particles show a broad distribution, from 
2 km/s (sensitivity barrier) up to about 25 
km/s, and averaging at about 10 km/s. These 
results are in agreement with results from 
the Pioneer 8/9 dust experiment (ref. 41), 
and the velocity distribution is expected as 
a consequence of the Poynting-Robertson ef- 
fect. 

In conclusion, one can state from the latest 
in-situ measurements that the Moon is hit 
by the apex-particle flux for particles > 1 /im 
in diameter at 10 km/s impact velocity and 
by the Sun-particle flux for particles < 1 /un 
in diameter at > 50 km/s impact velocity. 

Concerning the crater size distributions, 


Schneider at al. (ref. 6) have reported a 
bimodal distribution in crater sizes. Let us 
assume that the micron-sized and larger cra- 
ters have been produced by apex particles 
and the submicron-sized craters by sun par- 
ticles. It is then possible to convert the 
crater sizes into the respective particle sizes. 
For 10 km/s impact velocity, D P /d is be- 
tween 1.4 and 2 for crater diameters D P be- 
tween 1 /nm and 1 mm. For 50 km/s impact 
velocity and crater diameters <1 pm, D/d is 
4 for glass targets and 7 for metal targets as 
shown in figure 6. The conversion into fluxes 
is then accomplished applying formula (1) 
to the respective production crater distribu- 
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Figure 17. — Comparison of micrometeoroid fluxes 
derived from lunar microcrater measurements with 
those obtained in various satellite in-situ experi- 
ments. 

tion. Since the Moon is a spinning Earth 
satellite, the surface is hit by the dust 
streams only part of the time. The results, 
therefore, have to be multiplied by the factor 
of it , if unidirectional fluxes are concerned. 

The final flux results are given in figure 
17. The numbers of dust particles are cumu- 
latively plotted as a function of particle mass. 
The results obtained from lunar microcrater 
statistics are labeled “Moon.” This is the 
time-averaged flux on a lunar surface ex- 
posed directly to the two major dust streams 
while the Moon rotates and orbits the Sun 
with the Earth (ref. 16). 

In comparison with other flux results given 
in the flux diagram, it seems that the Moon 
flux is lower. However, all measurements with 
exception of the HEOS value are calculated 
under the assumption of a 2-7r-steradian de- 
tector exposure geometry and an isotropic 
flux direction. As discussed above, this as- 
sumption has to be revised drastically. For 
a two-stream particle flux distribution, the 
HEOS flux comes close to the Moon flux. This 
suggests that the present meteoroid flux is 
nearly the same as that averaged over more 
than 10 3 yr. 

Results contrary to this have been obtained 
based on solar flare track exposure age mea- 
surements of individual lunar micro-craters 
(ref. 42). The age distribution of micro- 



SOLAR FLARE TRACK EXPOSURE AGE 11000 yrl 

Figure 18. — Solar flare track exposure age measure- 
ments of single microcraters (Dp >20 pm) on 
sample 15205,51. 

craters found in this way suggests that a 
higher flux exists now compared with that 
existing at any time during the last 100 000 
yr (fig. 18). Further measurements are 
needed to explain this discrepancy. 

Erosion of Lunar Rocks 

Erosion of lunar rocks is dominantly con- 
trolled by micro-meteorite impact, which is a 
random process. Analytical model calcula- 
tions yield only statistical averages. How- 
ever, for most applications in lunar rock 
analysis it is of interest how a rock sur- 
face topography develops under continuous 
meteorite bombardment ; deviations from 
averages and the uniqueness of single-spot 
exposure histories have to be taken into ac- 
count. Therefore, a Monte Carlo computer 
model has been worked out (ref. 43). A 
given test area was randomly covered with 
craters of 17 size classes. 

The total simulated crater diameter inter- 
val extends from 152 /xm to 2.4 cm in spall di- 
ameter. This range contributes dominantly 
to the abrasion of rocks sampled during the 
Apollo missions. 

The probabilities for the occurrence of dif- 
ferent crater sizes were calculated on the 
basis of production state crater distributions 
actually observed on lunar sample surfaces. 

The test area itself was subdivided into 
roughly 27 000 area elements (cells), which 
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were characterized by their center coordi- ample, after 100 000 craters have been 

nates. The computer kept track of the ero- produced, about 5 percent of the surface will 

sion state, i.e., the depth level of each area have been impacted not less than 11 times, 

element. Some of the main results of this Figure 20 illustrates the maximum devia- 
computer model are as follows. tions of certain coherent reference areas 

Figure 19 gives the percentages of surface from the average erosion depth of the whole 

area affected (i.e., number of area elements surface, again as a function of total crater 

affected) as a function of total numbers of numbers or model time. Depending on the 

craters produced. Notice that already after size of the reference area there may easily 

about 8300 craters 50 percent of the total occur deviations up to factors of 10 in both 

area has been affected, whereas it takes an- directions, which means that small sample 

other 70 000 to 80 000 craters to completely chips, usually available for laboratory inves- 

cover the remaining 50 percent. With in- tigations, cannot be considered as being rep- 

creasing model time the probability for one resentative for whole sample surfaces, 

cell’s being affected repeatedly by impacts If one applies different crater production 
is increasing, too. The curves in the lower rates derived from model fluxes of the inter- 
part of the diagram represent these multi- planetary dust (figure 17), the model results 
impact percentages of the surface. For ex- can be related with absolute times, and aver- 

MODEL TIME 



0 20 40 60 80 100 0 20 40 60 80 100 

TOTAL NUMBER OF CRATERS PRODUCED I10 3 ) TOTAL NUMBER OF CRATERS PRODUCED (10 3 ) 


Figure 19 . — Percentages of surface area affected X Figure 20 . — Maximum deviation from the average 
many times (X = 1,2,3. ...) as a function of total erosion depth of the whole test surface plotted 

number of impacts produced. against model time for different reference areas. 


600 


COSMOCHEMISTRY OF THE MOON AND PLANETS 


age erosion rates can be determined. They 
range roughly between 0.07 mm and 0.7 mm 
per 10® yr. 

Erosion rates (single particle abrasion) 
can be calculated analytically. The mass ejec- 
tion M e per impact with mass m at 20 km/s 
impact velocity can be derived from a for- 
mula by Gault (ref. 44), resulting in M e = 
103.575 to m 33 Following a procedure by 
Gault et al. (ref. 44), erosion rates of 0.1 to 
0.3 mm per 10® yr are calculated (ref. 17). 
This value is dependent on the crater produc- 
tion rate used. In this case, a rate consistent 
with the data in figure 17 is used. 

Chemistry 

The glass lining of lunar microcraters has 
been investigated for admixture of projectile 
material with an electron microprobe. De- 
tailed analyses of several craters on the glass 
coating of sample 60015,35 and of one crater 
on the crystalline sample 12002,100 have been 
performed. 

The coating of 60015 was initially glass a 
few millimeters thick, but has devitrified to 
elongate plagioclase crystals with 2- to 50-/* 
dimensions and pyroxene-composition inter- 
stitial material. Point electron microprobe 
analyses for Fe, Mg, Ca, Al, and Si of glasses 
lining several 100-/*-sized pits show no indi- 
cation of material in the pit glasses that is 
outside the range in composition of the host 
sample phases. Either the projectile material 
composition is similar to the host material 
composition, or little or no projectile mate- 
rial remained in- the pit glass in these cases. 

Similar measurements that suggest some 
presence of projectile material or fractional 
vaporization during the impact process have 
been performed by Mehl (ref. 45). 

A 0.5 mm pit on sample 12002,100 showed 
numerous mounds associated with its glass 
coating. The mounds consist of mainly iron 
and nickel, with an iron-to-nickel ratio of 
20:1. This indicates that the mounds are the 
remnants of an iron-nickel micrometeoroid 
that formed this crater (ref. 16). Figure 21 
shows the distribution of iron and nickel in 



Figure 21. — Nickel-iron mounds in a 5O0-/im crater 
on sample 12002,100. Microprobe secondary elec- 
tron scanning picture with corresponding Fe-Ka 
and Ni-Ka X-ray scanning images (frame width 

*=*80 pm ). 
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the mounds and the surrounding crater glass 
lining. 

Conclusions 

The hmar microcrater phenomenology has 
been simulated in laboratory experiments in 
almost all aspects. Additionally, the labora- 
tory experiments allow the conversion of 
crater diameters into projectile diameters 
and masses at given impact velocities and 
projectile and target densities. These experi- 
mental results are applied to lunar micro- 
crater populations. 

Various measurements of microcrater pro- 
duction size frequency distributions exist in 
different size ranges. A compilation of these 
data yields the size frequency distribution 
for the diameter range between 1 /jn and 
1 mm, with an uncertainty of up to a factor 
of 5. 

The time development of crater population 
is qualitatively understood. The crater pop- 
ulations start as production populations, 
where essentially no craters have overlapped 
and destroyed each other. As time goes on, 
crater destruction occurs by subsequent cra- 
ter superposition. The populations pass 
through a transition state, where the crater 
distributions are gradually deviating from 
the production distributions. Finally, an 
equilibrium state is reached: on the statis- 
tical average, the number of craters in this 
state does not change. Theoretical descrip- 
tions of the equilibrium state satisfactorily 
agree with the observational data. 

Surface exposure ages have been de- 
termined from particle track measurements. 
The flux of interplanetary dust impacting the 
surface of the Moon is obtained by corre- 
lating areal crater densities with surface 
exposure ages. Such a correlation observed 
for large craters (100 < D P < 1000 p.) on sur- 
faces approaching equilibrium with respect 
to cratering may be a basis for a flux de- 
termination or, alternatively, an exposure 
time measurement if the flux (or the crater 
production rate) is known. Another inter- 
pretation of such a correlation is that since 


both solar flares and microcraters occur at 
the exposed lunar surface, they both suffer 
the same effects of erosion, and the correla- 
tion observed reflects only the relative “erod- 
ability” of different samples. 

The latest in-situ micrometeoroid measure- 
ments are important for the interpretation 
of the lunar crater frequencies: two major 
dust streams are hitting the Moon from dis- 
tinct spatial directions — apex and Sun; the 
correspondent impact velocities are about 10 
km/s for the apex dust stream and about 
50 km/s for the Sun dust stream. 

Taking into account the spin of the Moon 
relative to the Sun, it is possible to calculate 
a unidirectional flux valid for the Moon. This 
flux is fairly consistent with recent in-situ 
measurements, considered under the same 
conditions (Pioneer 8/9 and HEOS 2 re- 
sults). On the contrary, Hartung et al. (ref. 
42) have found an increased flux at present 
from individual microcrater glass linings 
dated by means of solar flare track measure- 
ments. Further measurements are necessary 
to resolve this discrepancy. 

Monte Carlo computer simulation and ana- 
lytical treatment of the erosion of lunar rocks 
by micrometeoroid impacts lead to consistent 
erosion rates of several millimeters per 10 7 
yr. 

Measurements of the chemical composi- 
tion of the glass linings of several craters 
> 100 /tin in diameter yielded the unam- 
biguous presence of projectile material only 
in one case, where the remnants of an iron- 
nickel micrometeorite have been detected. 
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